Astronomy & Astrophysics manuscript no. paper 


© ESO 2009 


May 28, 2009 





The kinematics of NGC1333-IRAS2A - a true Class protostar 



C. Brinch^'^, J. K. J0rgensen^, and M. R. Hogerheijde' 



o 
o 

>>: 



oo 
< 

o 



> 
in 

in 
o 
o\ 
o 



Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
e-mail: cbrinch@astro.uni-borLn.de 

Argelander-Institut fiir Astronomie, Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany 



ABSTRACT 



Context. Low-mass star formation is described by gravitational collapse of dense cores of gas and dust. At some point during the collapse, a 
disk is formed around the protostar and the disk will spin up and grow in size as the core contracts because of angular momentum conservation. 
The question is how early the disk formation process occurs. 

Aims. In this paper we aim to characterize the kinematical state of a deeply embedded. Class young stellar object, NGC1333-IRAS2A, based 
on high angular resolution (< 1" « 200 AU) interferometric observations of HCN and H'^CN J = 4-3 from the Submillimeter Array, and test 
whether a circumstellar disk can be detected based on gas kinematic features. 

Methods. We adopt a physical model which has been shown to describe the object well and obtain a fit of a parameterized model of the 
velocity field, using a two-dimensional axis-symmetric radiation transfer code. The parameterization and fit to the high angular resolution data 
characterize the central dynamical mass and the ratio of infall velocity to rotation velocity. 

Results. We find a large amount of infall and very little rotation on all scales. The central object has a relatively low mass of 0.25 Mq. Our best 
fit is consistent with both the interferometer data and single-dish observations of the same molecule. 

Conclusions. As an object with a low stellar mass compared to the envelope mass, we conclude that NGC1333-IRAS2A is consistent with 
the suggestion that, as a Class object, it represents the earliest stages of star formation. The large amount of infall relative to rotation also 
suggests that this is a young object. We do however find the need of a central compact component on scales of a few hundred AU based on 
the continuum data, which suggests that disk formation happens shortly after the initial gravitational collapse. The data do not reveal a distinct 
velocity field for this 0. 1 Mq component. 
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1. Introduction 



Radiative transfer - circumstellar matter - ISM: kinematics and dynamics - ISM: individual objects: NGC1333- 



An open question in the subject of low-mass star formation is 
how and when the velocity field of the molecular gas in the 
protostellar envelope changes from being dominated by radial 
infall to being dominated by Keplerian rotation. In the early 
stages of star formatio n, gas falls toward the center as gravita- 
tional collapse sets in (I Shu 1977 ). Cores, however, are known 
to car ry a certain amount of rotational energy (Goodma n et al. 
19931) . so as a core contracts, the net angular momentum spins 
it up, flattening the infalling mate rial in the direction of the 
rotation axis ( Cassen & Moosman| ri981 : Terebev et al. 1984 : 
Adams et al E988tlBasull998l) . After the main accretion phase, 
the envelope is gone and a planet forming disk with gas in 
Keplerian orbits is left. 

In this paper, we study the kinematics of the deeply em- 
bedded NGC1333-IRAS2A protostar using subarcsecond res- 
olution data from the Submillimeter AiTay (SMA) as well as 
detailed two-dimensional line radiative transfer and the model 



framework presented by B rinch et al.l (120081) . In that paper it 
was shown that the evolution of the velocity field of a hydrody- 



namical simulation of star formation, from the collapse phase 
to the disk phase, can be approximated by a single parameter, 
namely the radially averaged ratio of infall to rotation. Under 
the assumption that the simulation used in that work provides 
a good approximate description of the formation of stars, this 
parameter may be linked to the age of a protostar and can thus 
tell us about the evolutionary stage of the system. 

Observationally, yo ung stellar objec t s are classified using 
the scheme laid out by Lada & Wilkin g (119841) with the addi- 
tion of the Class by Andre et al.l ( Il993h . In this scheme, the 
evolution of the young star is related to the shape of the spec- 
tral energy distribution with the addition that Class objects 
are defined as protostars with < Menv, which is equivalent 
to the original empirical definition of Class objects as stars 
with a high submillimeter flux relative to their bolometric lu- 
minosity. The mass of the central object is not easily measured, 
but one way to estimate it is by determining the velocity field, 
since the magnitude of the velocity vectors are determined by 
the central mass when the material becomes dominated by 
rotation. 
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NGC13 33-IRAS2 was first de tected through IRAS obser- 
vations by Ijennings et alJ (Il987h . Submillimeter continuum 
imaging has subsequently shown that it likely consists of 
three different objects, NGC1333-IRAS2A, -IRAS2B, and - 
IRAS2C, of which the two former are well isolated and de- 
tected at m id-IR waveleng t hs, e.g ., with the Spitze r Space 



Telescope (lJ0rgensen et al.l l2007bt iGutermuth et al.l l2008i) . 



These sources also sho w compact emissio r i in millimeter 
interfe rometric images jLoonev et al. 120001: |j0rgensen et al 



2004bl) . The third source, IRAS2C, shows relatively strong 
N2H^ emission po ssibly reflecting the a bsence of a cen- 



tral heating source (lJ0rgensen et al.l 12004b). Modeling of the 
high angular resolution millimeter contin uum observations 
dLoonev et alj|2003l : lJ0rgensen et al. I l2004bh revealed the exis- 
tence of a central compact continuum component in NGC 1333- 
IRAS2A, likely reflecting the presence of a central circum- 
stellar disk with a mass of ~ 0.1 Mq (depending on the as- 
sumed dust temperature and opacity). NGC1333-IRAS2A was 
targeted as part of the PROS AC survey o f Class embedded 
protostars ( J0rgensen et al. 2005bl 2007al). Th e 1" resolution 
submillimeter continuum data (J 0rgensen et al ]l2005b') showed 
that the central compact component was resolved with a diam- 
eter of approximately 300 AU. 

In this paper we investigate the kinematical properties of 
IRAS2A using a global parameterization of the velocity field. 
We are particularly interested in the velocity structure on scales 
<100 AU, traced by the high angular resolution Submillimeter 
Arra}0 data, where a circumstellar disk is thought to be present. 
The SMA instrument is ideal for this task because we can probe 
the warm, dense gas with high excitation lines at high angular 
resolution and thereby zoom in on the central parts of the ob- 
ject. We furthermore benefit from the fact that the interferome- 
ter is less sensitive to ambient cloud material. 

This paper is laid out as follows: In Sects. |2] and [3] we 
present the data that we use, and our model and best fit are 
described in Sect. |4] Discussion and a summary are given in 
Sects. |5]and|6] respectively. 

2. Observations 

Supplementing t he SMA compact con figuration data from 
2004 October 17 (ij0rgensen et al.ll200 5b'). NGC1333-IRAS2A 
was observed in the extended configuration on 2006 January 
8 with seven of the eight antennas in the array at the time. 
The weather was good and stable during the observations 
with a sky opacity at 225 GHz of about 0.05. Combined with 
the 2004 compact configuration data, the observations provide 
good (m, y)-coverage on baselines between 17 and 210 kA with 
a few baselines at 232 to 264 kA. The complex gains were cal- 
ibrated by observations of the strong quasars 3c84 (2.2 Jy) and 
3clll (2.8 Jy) every 15 minutes. 

The SMA correlator was configured with an identical setup 
to the compact configuration observations with 2 GHz band- 
width at 345 GHz (LSB) and 355 GHz (USB). For the fines 



' Tiie Submillimeter Array is a joint project between tlie 
Smitlisonian Astrophiysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded by the 
Smithsonian Institution and the Academia Sinica. 



predominantly discussed in this paper, HCN and H'^CN J- 
4-3 at respectively 354.5055 and 343.3398 GHz, the spectral 
resolution was 0.15 km s"' (HCN) and 0.30 km s"' (H'^CN). 
In the maps used in this analysis, we use the natural weighting 
scheme to optimize for signal-to-noise in the resulting image 
at a small cost of resolution. With natural weighting the RMS 
is 0.225 (HCN) and 0.147 (H'^CN) Jy beam"' in a 1 km s"' 
wide channel with a beam size of 1.2" x 0.8". The initial cal- 
ibration of the data was done with the SMA MIR package (lOi 



2007|) with additi onal image processing done using MIRIAD 



dSaultet al.lll995h 

In addition to these data, we use previously published HCN 
and H'^C N single-dish data from the James Clerk Maxwell 
Telescope (|j0rgensen et al.l l2004bl) as well as 3 mm interfer- 



ometric observations from the BIMA array ( |j0rgensen et al 
i2004ai) . 

3. Results 

Figure [1] shows the zeroth and first moments of the image 
cubes, integrated intensity and velocity, respectively for both 
the H'^'CN line and HCN line. In this figure, the contour lines 
show the zeroth moment and the grey scale shows the first mo- 
ment. The contours start at the three sigma level corresponding 
to 5.4 Jy beam ' for HCN and 2.3 Jy beam ' for H'^CN and 
increase linearly by three sigma in both panels. The emission 
is seen to be centrally peaked with only a bit of extended emis- 
sion at the 3-6 cr level in HCN. This is in good agreement with 
the (m, v)-amplitude plots in Fig. |2] where HCN is seen to be re- 
solved while H'^CN is more compact. A small gradient in the 
velocity field is seen in the north-south direction in the HCN 
moment map with a magnitude of a few km s ' over the 
1000 AU scales of the observed emission. A similar gradient 



was rep orted in emission on larger scales by |j0rgensenetal 
(l2004al) . The gradient is most likely associated with one of the 



two outflows driven by IRAS2A (see panel b) Fig.[T] — cor- 
responding to its direction and polarity. A similar gradient 
is not present in the H'^CN map. A fit of a linear gradi- 
ent to this map shows that the systematic velocity differ- 
ence across the source is at most 0.1 km s ' over a scale of 
2" and that the direction of the gradient fit is highly depen- 
dent on the clip level of the map and thus unconstrained by 
the data. 

When the two SMA tracks are combined we get a good 
(m, v')-coverage on baselines between 17 and 210 kA. This is 
equivalent to emission on scales of 0.9" to 12 .1" which, at the 
distance of NGC1333 (220 pc; Cernis 1990l) . corresponds to 
linear scales of 198 to 2662 AU. Some flux is also measured 
on scales between 232 and 264 k/i (167 to 190 AU), but at low 
signal-to-noise. Any emission on baselines shorter than 11 kA 
is filtered out by the interferometer Figure |2] shows the aver- 
aged (m, v)-amplitudes as a function of baseline between and 
210 k^. In panel a) is shown the H'^CN J = 4-3 fine. This 
panel shows that the emission on all scales up to 200 kA is 
largely compact with a strong increase in flux on scales below 
50 kA. The histogram indicates the expectation values if no sig- 
nal were present and the full line is our best fit model, which we 
discuss in Sect.|4] Panel b) has a similar layout, but for the HCN 



C. Brinch et al.: The kinematics of NGC1333-IRAS2A 



3 



4- Ul3 



W 'CN J=4-3 



0- 




-4 - 



I.O 



0.5 



0.0 



-0.5 



1-1.0 



a) H"CN J=4-3 



2 0-2 
Aa (arcsec) 



-4 



to 




Fig. 1. Zero (contours) and first (grey scale) moments of H'^CN 
and HCN J = 4-3. Contour lines showing the integrated inten- 
sity begin at 3cr and increase by 3cr=2.3 Jy km s"' (H^-'CN) and 
5.4 Jy km s"' (HCN), respectively. T he full Unes indicate th e 
directions of the outflows identified by Knee & Sandelil ( 2000b . 



J = 4-3 line. For this line the emission is also well resolved 
on baselines below 80 kA, but with a lower signal-to-noise 
on longer baseUnes, disappearing in the noise at baseUnes 
longer than about 160 kA. This is likely due to the HCN 4- 
3 emission becoming optically thick on small scales rather 
than HCN having an intrinsically dilferent spatial distribu- 
tion than H'^CN. 

Figure [3] shows three position-velocity (PV) diagrams of 
both H'^CN and HCN. The PV-diagrams show the emission 
contours along one spatial axis and a range in velocities. The 
three panels of Fig. [3] correspond to PV-diagrams along slices 
of different position angles. If only rotation is present, the PV- 
diagram along the axis of rotation should show no evidence 
of systematic motions, whereas the P V-diagram along the per- 
pendicular axis should. The three axes for which we show PV- 
diagrams are chosen so that one coincides with the direction 
of the velocity gradient seen in Fig. [T] (the same direction as 




100 150 
(uv)-distance (kA.) 



< 2- 





b) HCN J=4-3 


*\ 




- 








50 


100 150 200 




(uv)-distance (kX) 



Fig. 2. Averaged (m, v)-amplitudes over the velocity range of 
+3 km s ' from the systemic velocity of H'^CN and HCN. The 
histogram is the expected amplitude if no signal was present. 
The solid line is our best fit model as discussed in Sect.|4] 



the north-south outflow), one perpendicular to this (coinciding 
with the east-west outflow), and one which is inclined with 45° 
with respect to the other two axes (coinciding with the direction 
toward IRAS2B and IRAS2C). 

The first thing to notice is that the three position-velocity 
diagrams extracted for each of HCN and H'^CN along the 
different spatial axes appear quite similar. There is no direc- 
tion that obviously shows a clear rotation pattern and with the 
peak intensities at different velocities located exactly at zero 
off'set, there is actually very little sign of the velocity gradient at 
all, except for small skew in the contours in the left-most panel 
of H'^'CN and emission at the 3-6 cr level at positive velocities 
in the center HCN panel. The black lines that are also plotted in 
Fig.|3]show Keplerian rotation profiles for a central mass of 0.3 
Mq (full line) and 1.0 (dashed line) in the first quadrant. 
The curves in the two lower quadrants show free-fall veloc- 
ity profiles for the corresponding central masses. Obviously, it 
is not possible to discriminate between rotation and infall by 
comparing the emission to these curves, because the emission 
is very compact, even for HCN. The emission contours are, 
however, more consistent with the 0.3 curves than with the 
1.0 Mq curves. The second thing to notice in the PV-diagrams 
is how fast the emission drops ofi" along the spatial axis while 
the emission at small off'sets shows very broad lines. This is 
true for both HCN and H'^CN. It is therefore likely that the 
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Fig. 3. PV-diagrams of the combined SMA data sets of H'^CN 
J =4-3 (top) and HCN J =4-3 (bottom). The three panels 
correspond to slices at different position angles with PA=10°, 
PA=100°, and PA=145°. In the leftmost panels are shown 
Keplerian rotation (first quadrant) and free-fall (third and fourth 
quadrants) profiles for M,=0.3 Mq (full line) and 1.0 Mq 
(dashed line). The contour lines are shown as increments of 3cr 
(0.9 Jy arcsec - and 1.9 Jy arcsec ^ for H'^CN and HCN, 
respectively) and have the same values in all three panels. 

velocity gradient which is seen in the moment maps is due to 
the outflow rather than rotation. After all, this gradient does lie 
in the same direction as the ou tflow. These PV-diagram s also 



resemble the ones presented by Hogerheijde et al. ( 19981) for a 



number of known outflow sources. 

4. Analysis 

We describe the observed emission with the model for the den- 



sity and temperature from |j0rgensen et al.l (l2002h . who derived 
this model from analysis of the continuum image of the source 
at 450 and 850 fim as well as at its broadband SED from 60 yum 
to 1.3 mm. The density is described by a power-law, 

n(r) = 1.5 X 10''cm"^(r/1000AU)"'■^ (1) 

while the temperature was calculate d using a ID con t inuum ra- 
diation transfer code as described by' j0rgensen et al. (2002, see 
also Fig. If of J0rgensen et al. 200 4a). The oute r radius of the 
model is 1 .2 x 10"* AU. |j0rgensen et al.1 (l2004bl) also modeled 
the chemical abundance profiles using the integrated intensity 
of the optically thin H'^CN single-dish lines. They found that 
to reproduce the line ratios, a model with a depletion zone is 
needed. The resulting HCN abundance profile reads. 



X(HCN) 



2 X 10-** 
2 X lO-** 



< 7 X 10^ cm-3 
> 7 X 10^ cm-3 



(2) 



assuming an isotopic abundance ratio of '^C/'^C = 70. We fur- 
thermore add a central compact component, in accordance to 
J0rgensen et al.l (12005 al) who found a need for an abundance in- 
crease near the center. We describe this component with a jump 



in the abundance within the radius above 90 K (~100 AU) to 
X(HCN) = 7 X lO-'l 

The new addition to the model is a parameterization of the 
velocity field, which allows for both infall and rotation. This 
breaks the spherical symmetry of the model, making it two- 
dimensional, although we keep the spherical description of the 
de nsity and temperatu re. We follow the approach presented 
by iBrinch et al.l (l2008h with a velocity field parameterized by 



two parameters. 



CM. / - -V^sinff 
cos a 



(3) 



with M,, the mass of the central object, and the angle a as 
free parameters. The velocity parameterization is such that 
a = corresponds to Keplerian rotation and no infall, whereas 
a - njl corresponds to free fall toward the center and no ro- 
tation. The use of this two-dimensional velocity field further 
introduces two free parameters, namely the inclination of the 
rotation axis with respect to the plane of the sky and the posi- 
tion angle of the rotation axis in the plane of the sky. We add 
a mean turbulent velocity field through a Doppler /7-parameter 
of 0.2 km s ' to the model. This value has been chosen to be 
large enough to get smooth line profiles, but small enough not 
to broaden the lines significantly. While this parameter is some- 
what degenerate with the central mass in the sense that it affects 
the line widths, it cannot account for the different line widths 
in the interferometric and single-dish data. Tlie FWHM of the 
H'^'CN lines are 2.1 km s ' in the single-dish beam and 4.7 
km s^ for the lines meassured by the interferometer. Only 
by letting the turbulent field vary across the model could we re- 
produce this effect. We have however chosen to use a constant 
turbulent field to keep the model simple, and a small change 
in its value does not affect our result significantly. For an in- 
depth discussion of turbulent line b roadening we refer to, e.g., 
Ward-Thompson & Bucklevl (l200lh . 

We us e the molecular excitation and radia tion transfer code 
RATRAN (iHogerheiide & van der T^l2000l) to determine the 
level populations of the HCN molecule on a computational 
domain consisting of multiple nested regular grids. The grid 
resolves scales from 6 AU to 1.2 x 10^ AU, which is about 
the highest dynamical range of scales we can resolve with 
RATRAN without loosing the gain of the accelerated lambda 
iteration (ALI) algorithm which makes the code feasible to run 
on a standard desktop computer The resulting image cubes are 
afterward either convolved with the appropriate beam profiles 
for comparison to the single-dish lines or Fourier transformed 
for direct comparison to the visibility tables in the BIMA and 
SMA data sets. 

4.1. Model fit 

We run a grid of models in the parameter ranges M, e {0.0, 1 .0) 
and a e {0, 7r/2) and we compare the single-dish spectra and 
PV-diagrams from the SMA data of H'^'CN only. A ;i'^-surface 
plot is shown in Fig. |4]for an inclination of 40° and a posi- 
tion angle of 100°. We add the x' value of the two single- 
dish transitions, calculated channel by channel on a velocity 
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Fig. 4. H^^CN surface in the parameters and a. The cross 
marks the minimum which lies at M,=0.25 and ci'=1.10 
with ax^ value of 1 .84. The black lines are lines of constant ac- 
cretion rate, with lines corresponding to 1, 3, 5, 7, and 9 xlO"'' 
Mq yr ' going from the lowest line and up. 



range from -4 to +4 km s~\ to thex^ value of the PV-diagram, 
where we compare pixels with a signal-to-noise better than 
three only. The;^f^ is seen to be minimized for parameter values 
of Mt - 0.25 M0 and a - 1.10, which we take to be our best 
fit model. We have chosen to adopt an inclination parameter 
value of 40°, because this parameter is not well-constrained 
by our fit. For any value below 40°, we get an almost con- 
stant value and only for values higher than 40° we see a 
systematic increase in thex^- The reason is that our best fit al- 
ready favors an infall dominated model (i.e. almost spherical 
velocity field) and therefore the velocity field is mostly inde- 
pendent of the inclination. As we go toward models with more 
rotation, the fit rapidly deteriorates with increasing inclination, 
which only supports the result that little rotation is present in 
IRAS2A. Similarly, the position angle is unconstrained for the 
same reasons, in full agreement with the PV-diagrams shown in 
Sect. |3] We have also plotted lines of constant mass accretion 
rate in Fig. |4] The accretion rate is calculated as 



, 2GM, 
M = 47rroPo 



sin(Q') 



(4) 



where po is the mass density at the distance ro = 23 AU The 
radius at which we evaluate the mass accretion rate has been 



chosen to be the same radius as was used by J0rge nsen et al 



( l2004ah for direct comparison. Our fit is seen to be consistent 
with a mass accretion rate of 5 x 10"^ M© yr"^ 

The first fit, shown in Fig.|5] is that of the single-dish lines 
H^^CN J =1-0, 3-2 and HCN J =1-0, 4-3. Being optically 
thin, the H'^CN lines are quite sensitive to the abundance pro- 
fil e and they constrain the radial variations well, as discussed 
by lj0rgensen et alJ (l2005al) . The 7=1-0 lines show distinct hy- 
perfine components. We model these components under the as- 
sumption of LTE between the hyperfine levels. The HCN lines 
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Fig. 5. Best model fit to the four single-dish transitions of HCN 
and H'^'CN. 



have a much more distinct profile and they are much more sen- 
sitive to the velocity parameters. However, the large optical 
thickness of HCN and the possible contamination of outflow 
emission in these lines make it difficult to assess the quality of 
our model (not the fit) on basis of these lines. This is the rea- 
son why we do not include the HCN lines in our;^'^ analysis. 
That said, our overall best fit model is also the model which 
gives the best fit to the HCN single-dish lines alone and it re- 
produces the (m, v)-amplitudes averaged over + 3 km s ' from 
tlie systemic velocity nicely. Only the primary component of 
the H'^CN 7=1-0 line is considered when calculating the;^f^, 
but our best fit model produces a reasonable fit to the hyperfine 
components. 

Figure |6] shows the best fit H'^CN model PV-diagram 
(panel a) as well as two other models for comparison, M»=1.0 
Mq and ff=l.l (panel b) andM,=0.25 Mq and q'=0.2 (panel c). 
The contour levels in these plots are set to the same values as 
the levels in Fig.[3]and so these plots are directly comparable to 
the center panel (PA=100°) of that figure. The best fit model is 
seen to reproduce the features of the data well in terms of line 
intensities and emission distribution, while the model in panel 
b has too weak lines and the model in panel c shows a skewed 
emission distribution. The asymmetry seen in panel c is typical 
for a rotating velocity field and this feature is not seen at all in 
the data. A direct comparison between the SMA data and our 
model is shown in Fig.|2] Panels a) and b) show the fit to a sin- 
gle spectrum toward the center position of IRAS2A for H^-'CN 
and HCN respectively. The fit to the H'^CN spectrum repro- 
duces the entire line profile well. From the residual of the fit to 
the PV-diagrams (panel c) it can be seen that it is not only the 
center position that fits well. In the case of HCN, however, the 
best fit model does not work as well. While the red side of the 
spectrum is relatively well reproduced, also spatially, the blue 
side suffers from a too strong peak and a too narrow line wing. 
The line wing emission is most likely the outflow, similar to the 
wing seen in HCN 7=4-3 in Fig.|5] The strong asymmetry, pre- 
dicted by our model, is not seen in the SMA spectrum, where 
the blue and the red sides peak at the same intensity. Because 
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Fig. 6. Model PV-diagrams of H^^CN / = 4-3 at P.A=10°. 
a) Best fit model, b) High central mass model with M» = 1.0 
Mq and a^l.l. c) Strongly rotating model with M,=0.25 and 
a-Q.2. All three panels are directly comparable to the center 
top panel of Fig. [3] The contours levels are the same in all 
three panels and identical to those in Fig.|3l 
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Fig. 7. Fit to the SMA data. Panel a) and b) show a single spec- 
trum toward the center while panel c) and d) show the residual 
of the model fit to the PV-diagrams. 

we do see the asymmetry in the single-dish lines, the missing 
blue peak is likely due the HCN line becoming optically thick 
on large scales where the interferometer filters out emission — 
and where our model does not necessarily apply, e.g., due 
to departures from the simple spherical power-law density 
profile or associated temperature profile. 

4.2. Outflow 

IRAS2A is known to drive two strong outflows, along axes 
that are almost perpendicular to each other. We do not model 
these outflows, and therefore we do not expect to reproduce 
features in the data that are associated with the outflows. As 
an optically thin line, we expect that H'^CN is much less af- 



Fig. 8. A HCN / = 1-0 spectrum toward the center of IRAS2A 
as observed by BIMA. Two models are overplotted; our best 
fit model in solid line style and a model with M,=1.0 in 
dashed line style. The truncation of the model line wings at 2.5 
km s"' is due to the finite pixel size of our grid. 

fected than HCN, and indeed only the single-dish HCN lines 
show broad wings with excess emission not accounted for by 
our model. Both the single-dish and the SMA HCN spectra are 
contaminated by outflow emission, and if we go to the BIMA 
observations of the low excitation line J - 1-0, it immediately 
becomes clear that care needs to be taken when modeling the 
HCN lines. Figure|8]shows the BIMA spectrum toward the cen- 
ter of the source. Two models are overplotted, the best fit with 
a solid line and a model with Mt = 1.0 M© using a dashed line 
style. The observed spectrum is extremely broad with wings 
extending all the way to ±6 km s Our model does not repro- 
duce these velocities, and indeed, even if we increase the mass 
significantly, the fine width only increases a little. 

5. Discussion 

Within the framewor k developed on bas is of the hydrody- 
namical simulations (Brinch et al. 2008), we determine the 
central mass of IRAS2A to be 0.25 M© with a Icr uncertainty 
of +0.1 Mq. This value is considerably smaller than the mass 
of the envelope of 1.7 Mq obtained from modeling of the 
dust continuum observations of IRAS2A and in this sense, 
IRAS2A is a true Class object with M, < Mgnv It should be 
noted that the implicit assumption for the velocity field in this 
work is that material is free falling — if it is not rotating. It is 
possible that the circumstellar material is not in free fall, either 
because free fall velocity has not yet been reached or because 
material is being supported magnetically or otherwise. If this 
is the case, the mass parameter underestimates the true mass 
of the star by up to a factor of four, depending on how far the 
system is from free fall (Brinch et al. 2008). 

Figure|4]shows that our best fit is consistent with a mass ac- 
cretion rate of 5x10"^ M© yr"\ If the mass accretion rate has 
been constant, the current central mass would have been ac- 
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ere ted in only 5x10^ yr. However, the accretion rate may have 
been lower in the past or maybe most of the material does not 
accrete directly onto the star The age of the star, which we 
derive from kinematical arguments, is consistent with previ- 
ously estimated time scales for IRAS2A. Time scales based on 



the outfl ow dynamics have been estimated by J0rgensen et al 



(l2004bl) and lBachifler et alJ (11998.) to be 4x10^ and (3-7)xl0^ 
yr, respectively. lJ0rgensen et al. found a time scale of 

1.7x10^ yr from an inside-out collapse model fitted to single- 
dish observations, with an accretion rate about an order of 
magnitud e lower than what we fin d here. In contrast to the 
model of .Torgen sen et al. (l2004hh . which is derived from 
observations probing the larger-scale non-collapsing enve- 
lope, our high angular-resolution data probe the core kine- 
matics on scales of a few hundred AU. The difference be- 
tween these results therefore also suggests the Umit to the 
usabiUty of the simple self-similar collapse models for these 
earliest stages of low-mass protostars. 

If the bolometric luminosity of IRAS2A of 20-25 Lq 
originates solely from the accretion through the relation L - 
GMfM/ r, the material needs to accrete at a radius of approxi- 
mately 16-20 solar radii (0.07-0.09 AU), which is about a fac- 
tor of 6-8 bigger than t he estimated stellar radius of young 
protostars ('Stahler"l988') and within the anticipated radius of 
a disk. Still, the earlier interferometric studies of the con- 
tinuum emission of IRAS2A show that the larger-scale in- 
falling envelope cannot explain the continuum emission on 
scales of a few x 100 AU, but that a n additional com- 
pact component has to be pres ent there dLoonev et alJ2003l: 
|j0rgensen et al.' 2004b', 2005b). If this compact component 
indeed reflects the presence of a circumstellar ~ 0.1 Mq 
disk, it must also have been built up by the accretion flow 
within the lifetime of the object. It therefore appears plausible 
that matter is not accreted directly onto the central star. 

The question is why we do not see the rotation if the 
compact continuum component is interpreted as a disk. 
One possible explanation of the absence of a rotational ve- 
locity component is that the disk Ues exactly in the plane of 
the sky and we are viewing it face-on. This however, does 
not agree with the morphology of the two outflows which 
both have components in the plane of the sky. Another pos- 
sibUity is that it is not a stable disk yet but rather a pUe up 
of low-angular momentum material very close to the star, 
which is stiU dominated by infall. Such a structure could 
represent a "pseudo-disk" found in the models of collapse 
of magnetized cores bv Galh & Shu (1993). By contraction, 
this material could stiU settle in a rotating disk on a "lo- 
cal dynamical time-scale" - essentially the free-fall time- 
scale. At a radius corresponding to the SMA beam (0.5" 
= 110 AU) this time-scale is a few xlO^ yr, comparable to 
the time-scale found above. Such low angular momentum 
material may therefore not have settled in a rotationally 
supported disk, yet. Finally, since IRAS2A drives two out- 
flows, it is reasonable to assume that IRAS2A is a close bi- 
nary. If this is the case, the compact dust component may 
lie inside the tidal truncation radius, preventing material 
from settUng into a Keplerian, rotationaUy supported disk 
jArtvmowicz & Lubow..l994.) . 



StiU, the main conclusion remains that IRAS2A appears 

to be an example of a very young protostar based on our best fit 
velocity model, with both the velocity field parameters point- 
ing toward a young source. Its velocity field is characterized 
by a high degree of infall relative to rotation, which in com- 
parison to the study of the vel ocity field in hyd rodynamical 
simulations of star formation by lBrinch et al.l (l2008i) suggests a 
young source. In that paper it was shown that the a parameter 
evolves smoothly with time from a high value (n/l) to a low 
value (0.0) as a core collapses and a disk is form ed. A study 
similar to this one was done bv lBrinch et alj (l2007b for another 
young star, L1489 IRS, classified as a Class I source. For this 
source, a disk was detected based on similar kinematical argu- 
ments as we use here, and for that source an a value of 0.26 
was found. Based on the value of the a parameter alone, we 
conclude that IRAS2A is a younger source than L1489 IRS, 
and this is entirely consistent with the general perception that 
Class objects are younger than Class I objects. 



6. Summary and outlook 

In this paper we have modeled the velocity field of the Class 
young stellar object NGC1333-IRAS2A. We have constrained 
our model primarily by the PV-diagram oftheH'^CN emission 
observed at high angular resolution with the SMA. Despite the 
angular resolution we find no evidence for Keplerian motions 
on any scales, even though a very dense, compact component, 
interpreted as a disk, is needed on scales of a fewxlOO AU to 
reproduce the (u, v)-amplitudes, seen in the continuum data. 

Repeating this approach for a sample of protostars would 
allow statistical comparison of the velocity fields of these stars 
and maybe allow a new evolutionary classification based on 
a and Mt/Me,„. This will only be possible when the Atacama 
Large Millimeter Array (ALMA) comes online, because 
ALMA will be able to do observations like the ones presented 
here in snapshot mode as opposed to a full night per source 
needed by the SMA. Furthermore, ALMA will allow us to 
zoom in on the innermost parts of protostellar objects and tell 
us what happens to the velocity field as circumstellar disks 
form. In this paper we did not model the outflow emission 
that is clearly seen in the BIMA HCN 7=1-0 fine and also 
noticeable in the HCN single-dish lines. By using the SMA to 
observe high-excitation lines, we filter out a lot of the emission 
on large scale where we expect the outflow to be dominating 
the velocity field. Still, we do not fit the HCN 7=4-3 line 
very well. This illustrates the need of good radiation transfer 
outflow models which can be used to investigate the interaction 
between accretion and outflows on very small scales. 
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